Electron density stratification in two-dimensional structures tuned by electric field 
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A new kinetic instability which results in formation of charge density waves is proposed. A spatial 
period of arising space-charge and field configuration is inversely proportional to electric field and 
can be tuned by applied voltage. The instability has no interpretation in the framework of traditional 
hydrodynamic approach, since it arises from modulation of an electron distribution function both 
in coordinate and energy spaces. The phenomenon can be observed in thin 2D nanostructures at 
relatively low electron density. 

PACS numbers: 73.50.Fq, 71.45.Lr, 72.30.+q, 73.21.Fg 



Recent progress in microelectronics is related to great 
success in controlled fabrication of low-dimensional semi- 
conductor systems. That is why transport properties 
of semiconductor nanostructures, both of classical and 
quantum nature, lately attract wide attention. In this 
paper we discuss a new type of purely classical instabil- 
ity, which can be observed in 2D nanostructures. The 
instability results in a formation of charge density waves 
(CDW). The main feature of the phenomenon is its ki- 
netic nature. In contrast to usual current and density 
instabilities in semiconductors ||, the electron dy- 
namics in the kinetic instability can not be described on 
the basis of the local hydrodynamic parameters, such as 
electronic density, drift velocity, and temperature. Re- 
markably, the kinetic instability can develop even in the 
Ohmic regime, when the stationary value of electric cur- 
rent is proportional to applied voltage. 

Closely related phenomena of the striations formation 
are widely known in the gas discharge physics. The stri- 
ated discharge has been observed since M. Faraday and 
is regarded as one of the most typical discharge forms 
[@7 Hi H 0; H • In spite of this, a consistent theory of stri- 
ation is up to now absent. In the last decades it was real- 
ized that for description of typical striated discharges the 
fluid approach fails. It was shown that the hydrodynamic 
description is valid only for very high electron densities 
when the collisions between electrons are frequent enough 
for the maxwellization of electron distribution function 
(EDF) H, |l^. At lower electron densities, occurring in 
typical gas discharges, the phenomenon is of essentially 
kinetic nature. In this case, the EDF perturbation in the 
striations is varying both in space, and along the energy 
axis m |l3[ [l|, and it is impossible to pa- 

rameterize it in terms of perturbations of electron density 
and temperature. The kinetic striations mechanism was 
analyzed first in |lll n2, 13 . In it was argued, that 



the necessary conditions for kinetic stratification are: a) 
The momentum relaxation is much faster than the en- 
ergy relaxation; b) The energy relaxation is mostly con- 
trolled by the energy gain in the external field Fq and 
strong inelastic collisions with a large fixed energy trans- 
fer Wq ; c) There should exist a mechanism of a weak con- 



tinuous energy loss. In a spatially modulated potential 
U{z) = -Fqz + SU(z)^U{z) = SU{z + L), these condi- 
tions provide |l2|, |l^ for the resonant EDF response 
at i = Lo/m, where Lq — Wq/Fq, m — 1,2,... This 
"resonant" behavior corresponds to the widely known 
empirical Novak's rule ||, |, |, [l7|. An idea was put 
forward that under the conditions a)-c) the insta- 
bility develops which results in the formation of CDW 
with the same periods Lo/m. In |l2|, |l^, the elec- 
tron kinetics was analyzed only in a given fixed electric 
potential profile U{z). However, a complete analysis of 
instability requires self-consistent calculation of the po- 
tential perturbation dU{z,t) in terms of the carrier den- 
sities perturbations. Since the discharge field depends 
crucially on the ion motion and on the complex ion gen- 
eration processes, even a linear instability problem for 
the gas discharge plasma is still lacking a self-consistent 
solution. 

In this paper we will demonstrate that, in princi- 
ple, the kinetic stratification is also observable in low- 
dimensional semiconductor structures |p^ . Moreover, it 
turns out that for the semiconductors a relatively sim- 
ple self-consistent analytical solution can be found. The 
main simplification follows from the fact, that, in contrast 
to the gas discharge, a compensating positive charge is 
fixed and homogeneous. The stratification conditions a), 
b), and c) can be easily achieved in semiconductors. The 
momentum relaxation is usually fast compared to energy 
relaxation. The requirements b) and c) are also usually 
satisfied, the scattering by optical phonons with energy 
Wq and scattering by acoustic phonons working as strong 
inelastic and weak quasielastic energy relaxation mecha- 
nisms. We will show that effect can be observed in 2D 
quantum wells with small thickness. Spatial periods of 
arising CDW equal to Lo/m and can be tuned by ap- 
plied voltage. We assume that lattice temperature Tq, as 
well as Fermi energy, are small compared to Wq (in what 
follows for simplicity we put To = 0). The condition b) 
requires that electrons be " hot" , and their energies be of 
the order of Wq, i.e. the electron gas in this case is non- 
degenerate. We also assume that electron concentration 
is small and neglect electron-electron collisions. 
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FIG. 1: (a) Motion of an electron in homogeneous field, 
Ufiiz') = — Fqz. The difi'usive "staircase" trajectories slowly 
drift in z direction with velocity n/ Fq. (b) Motion in the mod- 
ulated potential U{z,t) — —Fqz + 5U{z,t). 



Let us consider the motion of an electron in an exter- 
nal field Fq assuming for a moment that the only scat- 
tering mechanism is elastic scattering. This leads to a 
diffusion in coordinate space. As far as energy relax- 
ation processes are " turned off' , the electron is infinitely 
heated by the field Fq diffusing over the kinetic energy 
W. Evidently, this diffusion is strictly correlated with the 
diffusion in the coordinate space, since the full electron 
energy E — W — Fqz is conserved. In fact, an elec- 
tron diffuses in {z, W) space along the fine E = const. 
Now, if we take into account sufficiently intensive optical 
phonon emission, the electron motion will be restricted 
by a shell < W < Wq. In the process of diffusion with 
a constant total energy E, the kinetic energy increases. 
Reaching the point W — Wq, electron loses the energy 
Wo and starts a diffusive motion with a lower total en- 
ergy E — Wq. a trajectory of the electron in the space 
{E, z) is shown on Fig. la. The interaction with acous- 
tic phonons plays a role of a friction force leading to 
the continuous loss of the electron energy with the rate 
K due to spontaneous emission (since we have assumed 
that Tq — 0). Later we will show that for a 2D system 
the rate of energy loss k does not depend on the kinetic 
energy. There are three different time scales in our prob- 
lem: a transport scattering time t, a characteristic time 



of electron heating by the electric field tq ^ Lq/Dq (here 
Lo = Wo/Fo, Do = Wot/M and M is the electron effec- 
tive mass), and a time Wq/k which characterizes a rate 
of energy loss due to the emission of acoustic phonons. 
We will assume that 



r < To < 



Wo 



(1) 



Inequality (^ provides that the acoustic phonon scatter- 
ing may be considered as a small perturbation. Due to 
this scattering the "staircase" diffusive trajectories move 
slowly down (along axis E) with velocity k. This means at 
the same time that the trajectories slowly drift with the 
velocity sq = k/Fq along axis z (see Fig. la). Since sq 
is inversely proportional to the applied field one can say 
that the motion of trajectories demonstrates a negative 
differential mobility. It is well known that the negative 
differential mobility should lead to the current instabil- 
ity ^. However, our case is more complicated than 
the usual Gann instability, because one has to follow the 
motion of diffusive trajectories instead of the motion of 
individual electrons. As for latter ones, their average 
drift velocity obeys the usual Ohm's law v — Fot/M and 
is much larger than so (since the inequalities @) may be 
rewritten as sq <C ^ y/Wo/M). The fact that the in- 
stability can be observed in the Ohmic regime indicates 
that the effect is purely kinetic and can not be described 
in terms of hydrodynamic parameters. 

As far as the elastic collisions are dominant (see 
Eq.(|l])), the EDF is almost isotropic (l9|, f{z, W, ip, t) w 
fi{z, W, t) + fa{z, W, t) cos(((5). Here fi is an isotropic part 
of EDF, fa cos{(p) is a small anisotropic correction, and 
Lf is an angle between electron velocity and applied field. 
Denote J{z,W,t) = yJW/AMfa- The equations for 
and J can be written as follows [|], |[ ^ 
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(3) 



where D(W) = Wt jM is an energy dependent diffusion 
coefficient. For simplicity we assume that r is energy 
independent and, consequently, Z)(VF) is proportional to 
W. The boundary conditions for Eqs. (||), (|^) read 

ji\w=Wo = 0, {FJ - K/i)|vK=o = F.J\w=Wf,- (4) 

Here F{z,t) = —dU{z,t)/dz, and U{z,t) is a poten- 
tial energy, which includes both the self-consistent po- 
tential created by electrons and the external potential 
Ua{z) — —Foz. Condition fi\w=Wo = corresponds to 
the limit of a very strong interaction with optical phonons 
(" black wall" condition) . The second boundary condition 
is related to the conservation of the number of particles 



3 



in inelastic collisions |20). Eqs. (^), @, and (^) at k = 
have a homogeneous stationary solution 



now _ no fWo 

•J — -Jo — -rrr- , Ji — JiO — ttt- in — - 
Wo Wo \W 



(5) 



where no is the stationary electron concentration (we 
assume the following normalization ftodW = no). 

Since JodW = nov, solution (H) corresponds to the 
Ohmic regime. According to we rewrite Eqs. (|^), 
(^ in variables {E, z, t), where E = W + U{z, t) is a full 
energy of a particle. The result is given by 



J^-D{E-U)^, 
oz 



dt dz ^'^ dt'dE' 



(6) 



(7) 



The motion of a particle in the space of new variables 
is restricted by the curves E = Ei{z,t) = U{z,t), E — 
E2{z,t) = Wo + U[z,t) (See Fig. lb). The boundary 
conditions (ffl) can be rewritten as 



fi\E=E2(z,t) — 0, 
(FJ - Kfi)\E=Ei{z,t) = FJ\E=E2iz,t)- 



(8) 



Since the total energy of an electron changes slowly (with 
the characteristic time Wo/k), it will be useful to intro- 
duce an electron density distribution over the E axis 



NiE,t) 



Z2{E,t) 



dzf.,. 



zi(E,t) 



(9) 



where zi(E^t), Z2{E,t) are the inverse functions of 
Ei{z,t), E2{z,t), correspondingly, and the value NdE 
represents the number of electrons on "staircase" trajec- 
tories restricted by E and E + dE (see Fig. lb). The 
stationary value of N is given by iVo = no/Fg. Introduc- 
ing the notation lE{t) = J(z2, E^ t) (the stationary of Ie 
being equal to Jo ) we find from Eq. (Q) 



J{zi,E,t) = lE{t) + dz 

J zi 



dt at ' dE 



(10) 



Taking into account that dzi/dE — —1/F{zi,t), 
dzi/dt = {dU/dt)/F{zi,t) and using Eqs. (|), (0) we 
obtain the continuity-like equation that governs the elec- 
tron motion over the axis of total energy 
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JdUV 
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= 1 



E+Wo 



{t)~lE{t). 



Here the angle brackets mean averaging over z 
'dU\ 1 dU 



dt 



N 



dz^fr{z,E). 



(11) 



(12) 



Next we consider the deviations from the stationary 
solution in the linear approximation. A small periodic 
over coordinate modulation of the potential U — Uo ^ 
SUq exp{—iujt + iqz) induces the energy dependence of 
quantities Ie, N in forms Ie — Jo — SIqeKp(—iLut — 
iqE/Fo), N - No = SNg c^j>{-iujt - iqE/Fo). We will 
demonstrate that for q ~ qm — ±27rm/Lo (where m — 
1,2, ..) the imaginary part of lj is positive which implies 
that a stationary solution (||) is unstable. For q = qm 
the solution is periodic function of energy with a period 
Wo /m and Ie+Wo (0 = (t) . Then linearization of Eq. 
(0) yields 



(13) 



Fo + AF„ 



where AF™ = —iq^No (SUm) /5Nm- We see that the 
physics of the problem is governed by the only parameter 
(SUm) /SNjn (the subscript m implies that all quantities 
are taken at q = qm)- This parameter has a transparent 
physical meaning of a response of the averaged potential 
with respect to a small variation of electron density in 
energy space SN^- The instability (Im(ajm) > 0) occurs, 
when Re (SUm) /SNm > 0. In order to find this param- 
eter one should go beyond the averaged kinetic equa- 
tion ( prj ) and solve Eqs. (^, (^) together with Poisson 
equation. As long as Eq. (03) is already proportional 
to a small parameter k, one can simplify the solution 
of Eqs. (|6|), (0) assuming that k = 0, and neglecting 
dfi/dt and dU/dt (since Eq. ( p^ ) provides that w ~ k). 
Then Eq. (|^) reduces to dJ/dz = 0, which implies that 
J(z, E, t) = lE{t). As a result, Eq. (i) yields 



Mz,E,t)^lE{t) 



dz' 



D{E - U{z',t))' 



(14) 



The small variation of the distribution function 5fi can 
be found by linearization of this equation with respect to 
dim, SUm, the hmctions zi and Z2 being also linearized. 
The Poisson equation gives us a proportionality between 
6Um and the Fourier transform Snm of the variation of 
electron concentration 



dn{z, t) ~ 6 



dEU 



(15) 



El 



Here the variation includes the variation of 5fi as well as 
variation of the integration limits Ei(z, t) and E2{z, t). In 
this paper we restrict ourselves to the case of 2D semicon- 
ductor quantum well, assuming that the dielectric con- 
stant e is the same both inside and outside the quantum 
well. For such structure 
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5nr, 



(16) 



Using Eqs. (|T^), ( [iq ) and linearized Eqs. (g), (14) one 
can find the relation between 5Nm and 5Um- To calcu- 
late the parameter AFm, one should also average the 
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FIG. 2; Instability increment as a function of q. Instability 
regions correspond to g ~ = ±2nm/Lo 



variation of potential 5Um exp(i(7mz) with the stationary 
distribution function fio{E,z) (since we solve the prob- 
lem in a linear approximation). After cumbersome but 
rather straightforward calculations finally we get 



w„i = Soq„i + i 



(y.'n 



Wol + \ma„ 



(17) 



where a™ = J^" dy{l-exp{iqmy))/y, A„ = e^no/meFo. 
It is easy to check that for any m, Im(a'm) > 0. Thus, for 
Q = <lm stationary solution is unstable. For a weak field, 
Am ^ Ij the increment is field independent, lin{ujm) ~ 
k/Wq. One can show that for q w g,„ the spectrum reads 



D* 

uj{q) = Wm + (g - q,n)v - 



qm)\ (18) 



where D* = Do 1 



An 



. , (we neglected 

small corrections of the order of k to w and D*\ This im- 
plies that instability exists only in a small vicinity of q^ 
(see Fig. 2). This instability should lead to the formation 
of CDW with the periods Lq/tti. 

Next we discuss a possibility of observation of the ef- 
fect. The instability increment is proportional to the rate 
of energy loss k, which can be calculated for electrons in 
2D quantum well in full analogy with the 3D case 
For the case of infinitely deep rectangular quantum well 
of width a, calculations yield 



pa^h 



(19) 



Here Cq is a deformation potential constant, p is density 
of the crystal. This result justifies our assumption that 
K does not depend on electron kinetic energy. Also we 
see that k rapidly increases with decreasing a. The law 
K ~ can be understood from simple estimates. The 
momentum transfer from electron to phonon in the direc- 
tion perpendicular to the quantum well is of the order of 
H/a. Emission of such a phonon leads to the energy loss 
~ hS/a, where S is the sound velocity. The energy loss 
rate by the emission of longitudinal phonons may be ne- 
glected due to a small factor k\\a, where k\\ is the in-plane 
wave vector of 2D electron. We find that k is proportional 
to the integral over dqz of the product of energy loss hS/ a 
by the squared matrix element V?' q 1/a. The upper 



limit of the integral is of the same order, of 1/a, yielding 
K ^ a~^. This implies that the instability is more likely 
to be observed in thin 2D structures. On the other hand, 
the instability is suppressed by the electron-electron colli- 
sions, which lead to maxwellization of the EDF. Thus, the 
instability condition is given by Im(ti;„i) > l/rge, where 
Tee is the characteristic time of the electron-electron scat- 
tering. Crude estimate of r^e for hot electrons with char- 
acteristic energy Wo gives ~ e^no/e^HWo- Having in 
mind Eqs. (|^), (|l^), one can see that for low electron 
densities e^riQ/e^h < k, a certain field interval exists, 
in which the instability can be observed. Simple esti- 
mates for GaAs and GaN show that for thin quantum 
wells, a « 30A°, the electron concentration is restricted 
by small but quite reasonable value ^ 10^°sm~^. 

In conclusion, we have presented a self-consistent the- 
ory of kinetic stratification. We have shown that the 
spatial periods of strata equal to Wa/Fam and can be 
tuned by applied voltage. 
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